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I. INTRODUCTION 


In the last decade, rapid advances have been made in the design and fabrication of 
two-dimensional materials^ which can be used to manipulate light at small scales.^ The 
highly active held of plasmonics focuses on the interaction of electromagnetic radiation at 
the mid- and near-infrared spectrum with the collective motion of electrons in conducting 
interfaces and nanostructures.^ A goal is to generate electromagnetic waves that propagate 
with relatively small energy loss close to the interface between a conducting material, e.g. 
graphene, and a dielectric.In plasmonics such lateral waves should decay fast enough away 
from the interface; while, on the other hand, they should attenuate slowly enough in the 
direction of propagation along the interface. A wave that has attracted much attention in 
this context is the surface plasmon-polariton,^’^^’^^ with a variety of reported applications^^’^^ 
including invisibility cloaking, photovoltaics^^ and nanolithography.^® 

Most recently, direct experimental evidence was provided for generating surface plasmons 
by placing a receiving resonant antenna near a graphene sheet.Motivated by this advance, 
our goal with this paper is to analytically study the generation of surface plasmons by 
current-carrying sources via a solvable model for a fundamental setting. To this end, we for¬ 
mulate and solve exactly boundary value problems for the time-harmonic Maxwell equations 
in the presence of vertical and horizontal electric (Hertzian) dipoles near an isotropic and ho¬ 
mogeneous conducting sheet between two isotropic and non-magnetic unbounded media. 

The underlying theme, wave propagation near boundaries, has been the subject of impor¬ 
tant studies for over a century;see particularly the systematic and extensive treatment 
of Ref. 26. In these works, approximation techniques are developed for radiowave propa¬ 
gation; these have offered valuable insights into the lateral electromagnetic waves traveling 
near the boundaries between media of very different indices of refraction.^® In the frequency 
band of plasmonics,^ however, additional considerations have emerged because of properties 
of novel two-dimensional materials used in microscale applications. For instance, a thin layer 
of graphene has a complex surface conductivity, with a positive or negative imaginary part 
depending on frequency and doping, and introduces a jump in the tangential component 
of the magnetic held across the interface. The analytical consequences of this discontinuity 
are largely unexplored.This view suggests that the associated lateral electromagnetic 
waves be studied in detail in the near- and mid-infrared spectrum. 
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In this paper we analytically address aspects of the following question. What is the 
structure of the waves generated by dipole sources on a thin conducting film? Our goal is 
to single out fundamental attributes of the field which are intimately related to the film by 
analyzing a minimal solvable model. Our tasks can be summarized as follows. 

• We explicitly represent all field components in terms of one-dimensional Fourier-Bessel 
(Sommerfeld-type) integrals, in the spirit of Ref. 26. Our derivations, focusing on the 
electromagnetic field itself, differ from the use of the (non-physical) Hertzian potential 
invoked, e.g. in Ref. 7, by which the field components are derived via successive 
differentiations. 

• By a generalized Schwinger-Feynman representation for a class of integrals,we com¬ 
pute all field components via fast convergent series of known functions such as the 
Fresnel integrals, when both the dipole and observation point lie on the plane of the 
film and the ambient media have equal permittivities. Our model is thus simplified, 
yet without obscuring the goal to analytically understand the role of the interface. 

• In accord with applications in plasmonics,® we further simplify the exact solutions 
when the surface resistivity (inverse of conductivity) of the thin layer is much larger in 
magnitude than the intrinsic impedance of the ambient space. Then, a few terms are 
retained in the series expansions for the fields yielding simple approximate formulas 
for all distances from the source. 

Our approach is based on systematically solving Maxwell’s equations in the spirit of 
Refs. 26 and 30. Thus, we avoid any a-priori plane-wave approximations. We recognize 
a specific type of lateral wave as intimately related to the surface plasmon-polariton of 
transverse-magnetic (TM) polarization of plane waves^’^^ via the contribution of a certain 
pole in the complex plane of the dual (Fourier) variable; see also the treatments of Refs. 7, 
8, and 27. For a horizontal dipole on a thin film in free space, when contributions related 
to TM polarization in principle may co-exist with contributions of transverse-electric (TE) 
polarization, our analysis reveals that the TM surface plasmon-polariton, when present, is 
accompanied by a wave expressed by Fresnel integrals.^® 

The analysis presented here, with focus on explicit, physically transparent expressions 
for the electromagnetic field in terms of known functions, differs in methodology from pre- 
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vious studies of waves in similar settings/’®’^^’^’^ For example, in Ref. 11, the authors review 
dispersion relations for plane waves in plasmonics for a variety of experimentally relevant 
geometries, without discussing effects of point sources. In Ref. 7, exact integrals are formu¬ 
lated for the Hertzian potential produced by dipoles in the presence of a graphene sheet; 
the electromagnetic held is then computed by numerical evaluation of integrals. In Ref. 8, 
a similar task is carried out more extensively, with numerical evidence that the surface 
plasmon-polariton of TM polarization, recognized as a discrete spectral contribution, may 
dominate wave propagation under certain conditions on the surface conductivity. In Ref. 27, 
the authors numerically describe the held produced by dipoles near a graphene sheet, dis¬ 
tinguishing a “core region”, where the electric held can be much larger than its values in 
free space, from an “outer region”, where the held approaches its values in free space. 

Our work expands previous numerical approaches^’®’^^ in the following sense. By focusing 
on a simple yet nontrivial model with a conducting thin him, we are able to derive closed- 
form expressions which explicitly separate the primary held of the dipole, produced in the 
absence of the layer, from the scattered held which is sensitive to the him conductivity, for 
all distances from the point source. This approach singles out analytic aspects of the wave 
produced by the point source that are intimately connected to the surface conductivity of 
the him, thus showing how the surface plasmon related to TM polarization can dominate 
propagation in cases of physical and practical interest. Our analysis lacks generality, since 
we restrict attention to the case where the source and the observation point both he on the 
plane of the thin him; nonetheless, we view our treatment as a step necessary for tackling 
the problem of radiation by a realistic current-carrying source (rather than an incident plane 
wave) placed on the material surface. The character of the wave produced by the source 
of course becomes important at short distances or high enough frequencies. Our work aims 
at illuminating the complicated structure of the held in a simple nontrivial setting. 

The present work illustrates analytic aspects of the tensor (or, dyadic) Green function 
for the geometry of a thin conducting him at hxed frequency. Naturally, the response to 
any imposed current-carrying source can then be derived by superposition. This task lies 
beyond our present scope. Furthermore, we do not pursue numerical computations of the 
helds. The explicit computation by asymptotic methods^® of held components when the 
observation point or the dipole is away from the layer is left for future work. 

The remainder of this paper is organized as follows. In Section II, we describe the bound- 
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ary value problem for Maxwell’s equations. Section III focuses on the derivation of Fourier- 
Bessel representations, known as the Hankel transform,^® for the electromagnetic held. In 
Section IV, we evaluate exactly the requisite integrals when the dipole and the observation 
point are on the him, which is placed in a homogeneous space. Section V discusses the held 
resulting from our computations for distances from the source that are comparable to the 
wavelength in free space. Section VI concludes our paper with an outline of open problems. 

Notation and terminology. The time dependence is assumed throughout, where 
u is the angular frequency. M is the set of reals, Z is the set of integers, and boldface 
symbols denote vectors in M®. We write / = 0{g) (/ = o{g)) to mean that \ f /g\ is bounded 
by a nonzero constant (approaches zero) in a prescribed limit. / g implies f - g = o{g) 
in a prescribed limit. Sfttc {^w) denotes the real (imaginary) part of complex w. The term 
“sheet” has a two-fold meaning as either a material thin him or, as a “Riemann sheet”, a 
particular branch of a multiple-value function of a complex variable; and the terms “top 
Riemann sheet” and “hrst Riemann sheet” are used interchangeably. We use the terms “TM 
polarization” and “TE polarization” in the context of waves produced by dipoles to indicate 
the presence of certain denominators, denoted by V and Q in the main text, respectively, 
in the Fourier representations of the corresponding helds; each denominator appears in the 
rehection coefficient for the TM- or TE-polarized plane wave incident upon the thin him. 
The terms “surface plasmon-polariton” ® and “surface plasmon” are used interchangeably. 

II. BOUNDARY VALUE PROBLEM 

In this section, we formulate the boundary value problem for Maxwell’s equations. The 
current density of the vertical unit electric dipole, shown in Fig. 1, is 

Jix,y,z) = = e^6ro , tq = (0, 0, a) , (1) 

where a is the distance of the dipole from the layer, a > 0, is the unit vector along the 
2 ;-axis, and S,. denotes the Dirac mass at point r. For the horizontal unit electric dipole, 
shown in Fig. 2, the current density reads 

J (x, y, z) = J^{x, y, z) = 5ro • (2) 

The him has inhnitesimal thickness and scalar surface conductivity a, which is in principle 
complex and cu-dependent.^ The him lies in the plane z = 0 which separates region 1, the 
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FIG. 1. Vertical unit electric dipole at distance a from planar thin conducting film. The infinitely 
thin film lies in the plane z = 0 , between region 1 (half space {z > 0 } with wave number ki) and 
region 2 {{z < 0} with wave number /C 2 ); and has surface conductivity a. 

upper half space {z > 0 } with wave number ki, from region 2 , the lower half space {z < 0 } 
with wave number /c 2 (Figs. 1 and 2). We assume that ^kj > 0 and ‘^kj > 0 (j = 1, 2), i.e., 
a lossy medium j, including the case with Qkj ^kj. Note that /c| = where /io is 

the magnetic permeability of free space, since the media are assumed non-magnetic, and ij 
is the complex permittivity of medium j; in practically appealing situations, this ij has a 
small imaginary part. 

The time-harmonic Maxwell equations for the held (Ej,Bj) in region j dictate that 

V X Ej = iuBj , (3a) 

V X Bj = —i{k‘^/u)'Ej + /ioJ • (3b) 

By Gauss’ law, V ■ B^ = 0 and V ■ E^ = 0 if {x, y, z) 7 ^ fq; however, these equations are not 
independent from (3) in the time-harmonic case and, therefore, are not ntilized here. Eqna- 
tions (3) are snpplemented with boundary conditions for the tangential components, viz.,^’^^ 

X (El — E 2 )|^^g = 0 , (4a) 

X (Bi - B 2 )|^^q = MoO-Eji , (4b) 

where E|j := {Ei-(e^-Ei)e^}|^^g = {E 2 -(e^-E 2 denotes the (continnous) tangential 

electric held at z = 0. Notably, condition (4b) expresses the physical property that the thin 
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FIG. 2. Horizontal unit electric dipole at distance a from planar thin conducting film. The film 
lies in the plane z = 0, separating region 1 from region 2; and has surface conductivity a. 

conducting film amounts to an effective surface current of density Js = aEn at z = 0] this 
is viewed as a free current density for the held outside the him. The boundary conditions 
for the normal components of (E, B) are redundant for the derivation of a solution; see 
Appendix A. In addition to (4), we impose the Sommerfeld radiation condition, viz.,^® 

= as r-)■ cx) (r = \/x^ + + z^) , (5) 

uniformly 'wvvjriiz^ 0, for each scalar component ^js (s = x, y, z) of the vector-valued 
held = Ej, Bj); r = x + y By + z B;, is the position vector in Cartesian coordinates. 

Equations (3)-(5) with (1) or (2) constitute the desired boundary value problem. 



III. FOURIER REPRESENTATION OF SOLUTION 

In this section, we derive one-dimensional integral representations for (Ej, Bj) for a verti¬ 
cal and a horizontal electric dipole. The starting point is the Fourier transform with respect 
to {x,y) of (3) and (4). Accordingly, let 


dj{x,y,z) = 


(27r) 




( 6 ) 
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where ^3 = Ej, Bj and is the Fourier transform of in {x, y), assuming that the integral 
converges in an appropriate sense. Consequently, by (3) the transformed variables obey 

^ g ^ ^ 

^V^jz ~ '^^3y ~ : 

^ g ^ ^ 

-|- -^—Ejx itoBjy , 

t^Ejy tTjEjx xujBj^ , 

and 

^ g ^ ^ ^ 

iyBjz -g^Bjy i(^hj/u)'^Ejx yoJx > 

g ^ ^ ^ ^ 

Bjx i^Bjz i{kj/oj^Ejy -|- yo-Jy , 

i^Bjy iTjBjx iikjjoj^Ej^ HqJz j 

where (Jx, Jy, Jz) = (0, 0, 6a) or ((5^, 0, 0) is the Fourier transform of or J^, respectively; 
6a = 6{z — a). Equations (7) are complemented with the transformation of boundary condi¬ 
tions (4) and radiation condition (5). 


A. Vertical dipole 


Consider Fig. 1. In this case, by symmetry we have Bj^ = 0 for j = 1, 2. The remaining 
held components can be expressed in terms of Bjy, as shown below. Equations (7) combined 
yield the differential equation^® 

(^-^ + l^fjBjy = i^lj,o6{z - a) 2 :eM\{ 0 }, 

where 

ft - (k] -e- nY'" . 3ft > 0 (j = 1, 2) . 

In compliance with radiation condition (5), we write 

+ ft>0), j) = C<e-‘'’« (2 < 0) , (8) 

where C'> and C'< are integration constants to be determined. In view of (7), the remaining 
held components are given in terms of Bjy by the relations^® 


R 

Djx — ’ 


( 9 ) 



(;^ ^ 0, a) . (10) 


9 ^ _ iuri dByy - uj + S 

k] dz ’ k]^ dz ’ A;2 e 

To determine C> and C< we resort to conditions (4), by which Ei^ = E 2 x, Eiy = E 2 y, 
Bix — B 2 x = fJ^oO'E 2 y, and Biy — B 2 y = —fioaE 2 x at z = 0. In fact, we need only apply 
the hrst and third (or, second and fourth) of these conditions since we have set Bj^ = 0 ab 
initio; the other two conditions are then satished. Thus, we obtain 

r' — 7? 

0 >- , 

gi/3ia 

C< = fiokl^ p ; 


where the factor TZm is dehned as 


7^ 


m 


kl/32 - kl^i - ujiJ,oa/3i/32 

V 


( 11 ) 


which is associated with the reflection of TM-polarized plane waves from the thin layer (the 
subscript “m” stands for “magnetic”, implying TM polarization); and the corresponding 
denominator, V, is^’^ 

B = ki/32 + • ( 12 ) 


The two-dimensional Fourier integrals for (E, B) then follow from (6) with (8)-(10). 

To reduce representation (6) to one-dimensional integrals, we resort to the cylindrical 
coordinates {p,(f),z) where x = pcos0 and y = psin0 (0 < 0 < 2ti), following Ref. 26. The 
cylindrical components of the held = Ej,Bj) include '^jp = jx cos (p +'^jy sin cj) and 

= —'^jx sin (p +'^jy cos (p. Accordingly, let (^,?7) e-)■ (A, 0') with (^, 77 ) = (A cos 0', A sin 0') 
where A > 0 and 0 < 0' < 27r; thus, + yy = Apcos(0 — 0') and d.^ dp = Ad0'dA. By 
direct integration in 0', we hnd 


Ej^ = 0 and Bjp = 0 , 

as expected by symmetry because of the dipole orientation. For the remaining components, 
we invoke the known formula^^ 

E^^Jn{w) = - / (13) 

271 - Jo 

for n = 0 and 1, where Jn is the Bessel function of nth order, and n G Z. 

Consequently, after some algebra, we hnd the following integral representations. 
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For z > 0 (region 1) with {x,y,z) ^ (0,0, a), 

/ CXD 

dA A^ Ji(Ap) [sgn (;2 - , (14a) 

Bi4p.'^.«) = -^ rd^^^(V)[e‘'’‘''-“'-K»e"’'<'+“>] , (14b) 

Bi*(p. = r ! (15) 

and for 2 : < 0 (region 2), 

dXX^MXp) ^ , (16a) 

E 2 .(p,<t>,z) = (°°<lAA3j„(Ap)i , (16b) 

B 24 ,(p, P, z) = dA AOi(Ap) 1 . (17) 

In the above, sgn( 2 :) = 1 if > 0 and sgn( 2 :) = —1 if 2 : < 0; and, by (11) and (12), TZm 
and V are functions of A with 

ft(A) = (P|-A2)'12 (i = l,2). (18a) 

In accord with the Sommerfeld radiation condition, the top (physical) Riemann sheet in the 
four-sheeted A-Riemann surface for the held components is hxed by imposition of 

^/3j{X) > 0 , (18b) 


for each j = 1, 2; thus, /3j{X) = /3j{—X) in this Riemann sheet. We note in passing that in 
the limit where a J, 0 and z ^ 0, integrals (14)-(17) approach expressions that are divergent 
in the conventional sense yet become meaningful as hnite in the sense of Abel.^^’^^ This 
physically transparent interpretation permeates Section IV. 


B. Horizontal dipole 


Next, we focus on the geometry of Fig. 2. In this case, all cylindrical held components 
are in principle nonzero. By (7), Ejx and Bj^ satisfy^® 

g^2 + j j^2 ’ 

N 

32 


— -f- 13j ) Bj^ = 0 , 
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where 2 : 7 ^ 0 and /3j(A) is defined in (18), with admissible solutions 

^ (^ > 0) > ^ 2 x = (;^ < 0 ) , (19) 

Bi. = (;^ > 0 ) , (;^ < 0 ) , ( 20 ) 


consistent with radiation condition (5). Note that C'> 7 ^ C'< here because, in view of 
condition (4b), the tangential component of the B field is not continuous across = 0, in 
contrast to the formulation of Ref. 26. The integration constants C'<, C'>, and are 
determined through boundary conditions (4). By transformed Maxwell equations (7), the 
remaining field components can be expressed in terms of Ejx and Bjx as 


1 


- p c2y ^ 


/c • ^ - 


jx I , Bj^ _ _ ^2 


/c? - ^dB,^ 

u oz 




. . ^ dBjx 
ir]^Ej^ + u-^ 


B - e 


dEix 
oz 


JX 


(z ^ 0 , a) 

7 s V / -J \ oz / 

We proceed to compute the fields. By (4a) we impose Ejx = E 2 X at z = O; thus, 

unoikf - ^ 2 ) 


K^ = K, 


2(5ikl 


( 21 ) 

( 22 ) 


We henceforth set K = for ease of notation. By (4b), we enforce Bjx — B 2 X = 
at 2 : = 0 which yields 


C> = 




C< . 


Now introduce C = C'<. Thus, it suffices to determine the integration constants K and C. 

By (4), the continuity condition Ejy = E 2 y and the jump condition Bjy — B 2 y = —^oaE 2 x 
at z = 0 lead to the system 





/ ufj,oaf3i\ 

V ki-e) 


ki - e 


K + u 


/Si 


kf-e 


1 +^]+ 


/S 2 


ki-ej ki-e. 


(7 = 0 , 


kl^j kl^2 ^ 

+ 7^5 - 77 + OJllQa 


k'i-e ki-e 


2c2 


1 + 


Ti 


{kf-em-e) 


UT]^ 


1 + 


u}^oa/32\ 


kj-ev ' ki-ej ki-e 

After some algebra, the determinant for this system is found to be 

VQ 


C = ufi^e 


K 


i(3ia 


V = 


(kf-em-e) ’ 
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where V is defined by (12), and 


Q — /3i + /32 + 


(23) 


corresponds to TE polarization/ It follows that 


K = —ojjjiQ 


/i(^2 - + l^2{kl - ^2) + u^ioal3il32 

VQ "" 


Accordingly, all field components are now obtained via (19)-(22) by use of Fourier inte¬ 
gral (6). By representing the resulting two-dimensional integrals in cylindrical coordinates 
{p,(f),z) and using (13),^® we obtain the following expressions. 

For z > 0 (region 1); 

^ip(p,0,^) =-^cos^f^ [ dAA|[Jo(Ap) - J 2 (Ap)] 


dTT 


ki 


-[Jo(Ap) +J2(Ap)]-7^eje*^^(^+“) 

+ - J 2 (Art] 


(24a) 


Ei^{p,(f),z) = 


^^sin0^^ dAA{-[Jo(Ap) - J 2 (Ap)] ^7^e 

+ [Jo(Ap) + J 2 (Ap)] A 7 ^^|e*^d^+a) 
+ ^“dA a|a;2Jo(Ap) - y [Jo(Ap) + J 2 (Ap)]}^ 


(24b) 


Fi^(p,0,^) = ^^cos0 j dAAVi(Ap)[sgn(^ - a)e*^d^ ay g*pi(2+a)j ^ ^^240) 


Bip{p,(l>,z) 


ho . 

-sm 

Stt 


^ dAA{[Jo(Ap) +J2(Ap)]7^„ 


-[Jo(Ap)-J2(Ap)]7^e}e*^d^+a) 

N 

/OO 

dAAJo(Ap) 


(25a) 
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oo 




/^o 

Stt 


cos (j) 



dAA{[Jo(Ap)-J2(Ap)]7^^ 


-[Jo(Ap) +J2(Ap)]7^e}e*^^(^+“) 

N 

/ OO 

dAAJo(Ap) 


(25b) 


fii,(p,0,;^) = -^sin^^^dA AVi(Ap) ; (25c) 

and for 2 : < 0 (region 2), 

E2p{p,^,z) = -^COS 0 ^ dAA|[Jo(Ap) + J2(Ap)] 

+ [JoM - J2iXp)] ^ (26a) 

^20(p,0,^) = ^sin0^ dAA|[Jo(Ap) - J2(Ap)] 

+ [Jo(Ap) + J2(Ap)] ^ ( 26 b) 

E 2 z{p,(t),z) = -^^COS0 [ dAAVi(Ap) ^ , (26c) 

ZT^ Jo E 


B2p{p,J),z) = ^Sin0^ dAA|[Jo(Ap) + J2(Ap)] 

+ [Jo(Ap) - MXp)] ^ (27a) 

B 2 ^{p,(j),z) = cos (f) dAA|[Jo(Ap) - J2(Ap)] 




+ IJoM + J 2 (Ap)] Ije*'’--*' , 

(27b) 

B2z(P, (/>,z) 

ipo . 

= -sm 

27r 

700 1 

0 / dAAVi(Ap). 

Jo Q 

(27c) 

In the above, the factor 

TZe = 

02 — 01 + ujfXQa 

(28) 


Q 


is associated with the reflection of plane waves with TE polarization from a thin layer (the 
subscript “e” stands for “electric”, implying TE polarization); also, recall (11), (12) and (23) 
for TZm, V and Q, respectively. Similarly to the vertical-dipole case (Section III A), the top 
Riemann sheet is dehned by Q/3j(A) > 0, as in (18b). 
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C. Poles in A-Riemann surface 

We now discuss the role of singularities in evaluating integrals (14)-(17) and (24)-(27). 
The singularities of the integrands include: branch points at A = ki, k 2 because of the 
multiple-valued /3j{X) of (18a); and simple poles identihed with zeros of V{\) and Q(A) 
by (12) and (23). The A-Riemann surface associated with each integrand consists of four 
Riemann sheets if ki ^ k 2 ] the physically relevant (top) one is hxed by ‘^I3j{\) > 0. The 
location of the poles in this Riemann surface is closely related to the complex surface con¬ 
ductivity, a 7 

1. Denominator V{\) 

By (12), the zeros of V{\) satisfy 

{kf - A2)V2 + (A;2 _ A2)l/2 + ^= 0 . (29) 

Equation (29) is identihed with the dispersion relation for the surface plasmon in the context 
of TM polarization for plane waves^ if A is viewed as the component of the vector wave 
number tangential to the interface. For nearly lossless ambient media {Qkj ~ 0), roots 
of (29) lie in the hrst Riemann sheet if 

SJa > 0 . (30) 

For a graphene layer, this is consistent with the Kubo formula in the far-infrared regime. 

Consider the special case with ki = k 2 ='■ k (Section IV) in which the A-Riemann surface 
for the helds has two sheets. Then, 'P(A) = {2k‘^ + ufioay/k‘^ — \‘^)\/k‘^ — A^ and (29) reduces 
to yJkX — = —2k7/{ujjjiQa)^ which is solved in the hrst Riemann sheet if —tt < arg/c^ — 

argcT < 0. Equation (29) has two solutions, A = i/cm, where 



the branch of the square root is chosen so that Xkm > 0. These i/cm are simple poles of 
the corresponding integrands for the helds. In the “nonretarded frequency regime”,® one 
imposes \u}^Qa\ \k\ by which l/|a| 3> where is the magnitude of the 

intrinsic impedance of the adjacent homogeneous space; thus, we obtain km ~ i2k'^/{ufXQa). 
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For arga = 7 r /2 — 5, 0 < 5 < tt, the analytically continued square-root function yields 
= ie'^^2k‘^/{ufj,o\a\) and km lies in the hrst Riemann sheet for the appropriate 
range of S and phase of k, as outlined above; this km reads 

2k^ 

km ~’ l^hoo-| < \k\ . (32) 

uno\a\ 

If ‘^^yk‘^ — km < 0, then km is in the second Riemann sheet. 

We adhere to the following dehnition for the purposes of the analysis in Section IV. 
Definition 1 (TM surface plasmon). The TM surface plasmon corresponds to the 
residue contribution to the electromagnetic field from the pole \ = km, provided this pole lies 
in the first Riemann sheet. 

We note in passing that in the more general setting with ki 7 ^ /c 2 , one hnds the simple 
poles i/Cm ±i{k‘l + kl)/{(jJiJLQa) if \{kl + A;|)/(a;/ioO')| > max{|A;i|, \ k 2 \}. 

2. Denominator Q(A) 

By (23), the zeros of Q(A) obey 

{ki — + {kl — + ujnoa = 0 , (33) 

which has two roots, i/Ce, in the A-Riemann surface, where'^’® 

For definiteness, the branch of the square root is chosen so that ke ^ kif one sets ki = k 2 = k 
and then lets oonoa —)■ 0. The i/Ce are simple poles of the corresponding integrands for the 
helds. By (33) and for nearly lossless ambient media {Qkj ~ 0), ±/Ce lie in the top Riemann 
sheet if 

9a < 0 , (34) 

in contrast to (30). Thus, for the appealing case with nonzero 9a and an ambient lossless 
dielectric,®’^ only one of km and k^ lies in the top Riemann sheet.^ If |a;/ro<^P ^ l^i + 
then ±ke ~ ±{i/2){k1 — k1)/{ujj,oO') for ki Y ^ 2 ] hence, \ke\ may take a wide range of values. 

In the special case with ki = k 2 = k (Section IV), (33) reduces to — A^ = —ufXQa/2, 
which is solved in the hrst Riemann sheet if —tt < arga < 0. We hnd 

= ( 35 ) 
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If |a;/ioO'| \k\ then kg ^ k — {uj{Sk) ~ k. By setting arga = 7r/2 — S, we find 

kg — k ^ ^ , \u:^Qa\ \k\ (0 < 6 < tt) . (36) 

8k 

The above discussion leads to the following remark. 

Remark 1 (On the co-existence of poles km and kg in top Riemann sheet). Consider the 
case with ki = k 2 = k. We henceforth assume that 0 < arg/c^ < tt/ 2 and Sfta > 0 with 
arga = n/2 — 5 (0 < 5 < n). Thus, at most one of the poles at km and kg is present in the 
top X-Riemann sheet. 

Specifically: If'^cr > 0 with 0 < 5 < 7r/2 — arg then km lies in the top Riemann sheet; 
and if< 0 (n/2 < S <n), then kg lies in the top Riemann sheet. 


IV. EXACT EVALUATION OF FIELD ON PLANE OF FILM 

In this section, we analytically evaluate the integrals for the electromagnetic held (Sec¬ 
tions III A and IIIB) on one side of the layer by setting ki = k 2 = k and allowing the dipole 
and observation point {x, y, z) to approach z = 0. 


A. Formalism 

For the vertical dipole (Fig. 1), let a 0 with 0 < 2 : < a, so that (p, 0, z) approaches the 
him from region 1. Equations (14) and (15) reduce to the following expressions. 


(vd) lUfio 

Eip - 


^(vd) _ -r-u 
1</) 


27r (ca/ioo-)^ [ p 
2 


2kf upoa 

- h — 


ipo / Jf_ 

2k^kl,p^ 


-T - —T 
2ti \ (jjpoa p ^ ^ 


P 

k^ 


2 I 2 k^p *^4(^mp) 

2k^ 


{ujpoaf 


j , k^kf^p 
I2 H- HykmP) 

UjpoCJ 


{upoaY 

(vd) _ ^l^(^R(vd) 

~k^pdp^P 


Ri^kmP^ 


upo Id f 1 1 

-- Jl- 7:^2 + 


2TikP p dp^\upQa p 


P 


k^ - 


2k^ 


k'^klp 

I2 H- IzykmP) 


{ujpoaY\ upoa 


‘^k^kl^p^ 

{upoaY 


RikmP) 


(37) 


(38) 


(39) 
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where the requisite integrals read 



(40a) 

(40b) 

(40c) 

(40d) 


with — {kpY]~^[l 2 + Thus, we only need to compute / 2 , / 2 , h and 14 ^. 

For a horizontal dipole (Fig. 2), we choose to evaluate the held on the side of the layer that 
faces the source-free region: Let a 0 and t 0 so that the observation point approaches 
the plane z = 0 from region 2. By (26) and (27), the ensuing representations for the held 
components at 2 : = 0 in region 2 read 


p(hd) 

^ 2 z 


iupo 1 
27r {upoay 


coscj) 




+ 2k‘^kl^ph{kmp) 



4-k^ o 1 1 

- f^h{k,^p) + 2k^h{k^p) , (42) 

ca/ioo- J J 

dr = ^ 

- ^ hikmP) + 2k^ hikrnP)] ] , (43) 

upoo- J J 

^(hd) ^ ^ ^ /2 + upoaklp h{kep) - 2kl h{kep) , (44) 

OTT ]_ P P 

^(hd) ^ ^ [Zk"^p I^[kmp) - up^a h{kmp)] + -ph 

_ (-45) 

<’ = £ 4[2PV/3(fc„p) -u.p„,/,(fc„,p)] + 

- ph{kep)+ujpo(T h{kep) j. (46) 
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In the above, we provide expressions for the helds on one side of the layer for each source; 
the helds on the other side can be found through the imposed boundary conditions (Section II 
and Appendix A). The integrals hi^w), and are evaluated in Section IVB. (The 

reader who wishes to skip technical derivations may proceed directly to Section IV C.) 

The held of the vertical dipole bears TM polarization, as indicated by (37)-(39) with (40). 
On the other hand, the held of the horizontal dipole carries the signatures of TM and TE 
polarization. 


B. Key integrals and generalized Schwinger-Feynman representation 


Next, we provide derivations for requisite integrals (40b)-(40d) in terms of known tran¬ 
scendental functions. This section forms a crucial part of our formal analysis. We follow the 
formalism of Ref. 30, elements of which we repeat here for the sake of completeness. This pro¬ 
cedure should clarify the role of certain singularities related to material properties of the him. 

For later algebraic convenience, set kp =: iq and w =\ i?> recalling that w = kmP or keP] 
and treat q and h as positive for the purpose of computing the integrals. The extension of 
the requisite integrals to the actual, physically relevant values of these parameters will be 
carried out once the integrals are evaluated. Accordingly, set a/ (kp)'^ — A'^ = 
and treat a/A'^ -|- q^ as positive if A' is real. Firstly, we directly hnd that^'^’^® 

^ 1 _ ikp 

h = i{q + e-'^) = kp + =-tq-\l - e-^) = — - . (47) 

kp 


1. Integral 


By (40c), we resort to writing Ji(A) = (l/2)[i7j^^(A) -|- where is the (hrst 

or second) Hankel function of order z/.^^ By the identity = h[^\X), valid for 

A 7 ^ 0, (40c) becomes 


, hP(X) 

A'2 + B2 


= 2 


- lim 

2 e4,0 


hP(X) _ T „m, 






(48) 


by contour integration, where the contours C and are depicted in Fig. 3; cf. Ref. 36. The 


integrals in the hrst line of (48) provide the Cauchy principal value of 4 dA' i7|^^(A') (A'^-|- 
B^)“^ which accounts for the singularity of i7|^^(A') at X' = 0. 
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FIG. 3. Branch-cut configuration and contours in the A'-plane for the integral see (48). The 

integrand has simple poles at ±i6 (fi = —iw is treated as positive). Contour C = (—00, —e) U U 
(e, 00) (solid line with arrows) is identified with the real axis except for a small indentation, of 
radius e around 0 in the upper half-plane; e | 0. Contour C+ (dashed line) in the upper half-plane 
closes the path; the contour C U C+ picks up the residue at A = ffi. 

2. Integrals Ia{w) and Ii{w) 

We seek exact formulas for (40d) with w = kmP and w = kep. Such formulas should be 
particularly useful for \upQa\ \k\, by which \ke\ ~ \k\ \km\'^ cf. (31) and (35). Thus, we 

view as small the parameters: \k/km\, with zu = kmP', and \ {ke — k)/ke\, with zu = k^p. Our 
derivation is tailored to the case where km lies in the hrst Riemann sheet (see Remark 1); the 
extension to other values of the physical parameters is discussed in the end of this section. 

The hrst step is to convert Ii{w) in (40d) to an integral of an elementary function. In 
the spirit of the Schwinger-Feynman approach, consider the generalized representation^® 

1 poo POO p1 

/ / dsi ds2 = u duu-^+^ [(2t-C)M+C]-^-‘^, 3fJz/ > 0 , 

r(z^) Jo Jo Jo 

via the change of variable (<§1,52) e-)■ {u,v) with (si,S2) = {uv, (1 — u)v) and integration in 
v] r(; 2 ) is the Gamma function. We assume that 0 < C < 21, a condition to be relaxed later; 
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here, 21 = A'^ + € = A'^ + By using the above representation for 21 we obtain 


, (A'2 + q^) 


2 \ — v 


L A'2+B2 


Ji(A') = 


1 — e 


i2'KV 




Co 


r(i + z/) 


X [q^u + B^(l - u)] 2 *' Ki+^(^{q‘^u + B^(l - u)Y^‘^) 


(49) 


The contour Cq serves the analytic continuation of the integral on the left-hand side to values 
of z/ with 3fJz/ < 0 (see Fig. 4). To derive (49), we interchanged the order of integration (in 
A' and u), and made use of a result from Ref. 35 to carry out the integration in A'; is the 
modified Bessel function of order z/, with Ki/ 2 {C) = a/tt/( 2(C) By the change of variable 

n !-)■ C with = q^u -|- B^(l — n), integration by parts via the identity — 6^)'^] = 

— B^)*^"^, and the subseqnent substitntion z/ = —1/2, we find 

,, VA'2 + g2 




TdC 


c 


B2 


6 


= _T(l-e“'')+» 


ZU^ 


{kpf 


zu 


d?7 


.-c 


7] 




Jkp/m 

which provides Ia{'uj) in terms of an integral of an elementary function. 
We proceed to simplify (50). For this pnrpose, dehne 

7] 


(50) 


/ dr/ 




vC 


(51) 


/e \/ i - T 

We seek appropriate expansions of this hF for ^ = 1 with: (i) i = k/km and zu = kmP where 
\l\ < 1; and (ii) i = kjke and zu = k^p where |1 — f'| < 1. \zu\ should be unrestricted in each 
case. 

(i) Case £ = k/km, ^ = kmP- By writing 


IF {k jkm, Ij kmp') — hF(0,1, kmp') hF(0, k/km, kmp') 


(52) 


and resorting to Ref. 31, we obtain the formnla 

1F(0,1; kmP) = 1-2 - Yi^kmP)] + y H^i\kmP) , (53) 

where IHI/y (W) is the Struve (Neumann) function of order v. In (53), the last term cor¬ 
responds to the contribution to I^iJvmP) from the residne at A' = kmP- For the remaining 
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FIG. 4. Contour Cq of integration for generalized integral on right-hand side of (49) pertaining 
to / 4 (cu); 0 < . The integrand has branch points at u = 0, —— 6^); only the cut 

emanating from 0 is shown here, indicating the passage of the integrand from the first Riemann 
sheet for ^ where > 0 if u > 0, to an adjacent Riemann sheet. 


terms, note the large-|C| expansion^ 


H4C)-U(C) = --'Ef;hi^G 

^ r(z/+2 - n) VC 


2n—v-\-l 


+ (^(ICI 


-2M+U-1 


) , largCI < TT . (54) 


Now consider the term W{0, k/km] k^p) in (52) by inspection of (51). By expanding 
T]{1 — at T] = 0 with |? 7 | < 1 and integrating term by term, we derive the geometri¬ 


cally convergent series 


W (0, k/ /Cm, kyi^p') ^ ^ 


n fi-e< 


klJ V c 


^2 

(=kp '^m 


^ < 1 , (55) 


where {b)i is Pochhammer’s symbol; {b)i = r(/ -|- b)/r{b). Formula (55) concludes the 
evaluation of W{k/km, 1; kmp) by (52). 

(ii) Case i = k/ke, w = k^p- By changing integration variable from 77 to u = 1 — ?7 and 
expanding 77(1 — 77 ^)“^/^ = (1 — v)v~^^‘^{2 — at u = 0, we recast (51) to 


W{k/ke, 1 ; kep) = 


1\ l + 2l{2kepy 




fl{{ke - k)p) , 


where 


fi{z)= / du 
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This fi{z) has generating function 

®(C; z) := f; (1 + + C)^) 


1=0 


where J^o is given in terms of the Fresnel integrals C{z) and S{z),^'^ viz., 


-ir 

J^o(z) = I dr , = C{z) — iS{z) 


'0 


V 2 


TFT 


= - e-‘^Z(z) i J^(z) := e‘ 


g-*r 


dr 




'TIT 


In particular, note the asymptotic expansion 


Z(z') = 


1 + ^ + oWt 


\/2'nz 

Thus, by use of the generating function 0 we write 


as z —)■ cx) 


(56a) 


(56b) 


/i(r) = ^/Vo(r)] , 

and derive the desired expansion for W: 




1=0 


^+1/2 Jl 


2 = P 


(57) 


This series exhibits geometric rate of convergence for \{ke — k)/ke\ < 1 and all p where 

Z2S {.^l^oW\fp 


p := {ke - k)p ~ e" 


Sk 


\ujpoa\ < \k\ , 


(58) 


with arga = 7r/2 — 5 (0 < 5 < tt). This p is analogous to Sommerfeld’s “numerical 
distance”;^® cf. (36). 

Thus, by (50), integral Ii{w) reads 

Ii{zo) = —(l — e*^^) + ^ W{kp/w, 1; zu) , (59) 

zu^ zu 

where W is provided by (52) with (53) and (55) for zu = kmP] and by (57) for zu = kep. 
By (40d), Ii{uu) equals 


hiuj) = 


^-^1(1- + i {kp )_ (^kpluu, 1; w) 


1 kp 
kp uu"^ 


uu 


(60) 


zu"^ — {kpY 

In the above, ‘R^Juu‘^ — (kp) 2 = ^^{kpf — uu"^ > 0, in accord with the dehnition of the top 
A-Riemann sheet, and km is assumed to lie in this Riemann sheet; see (32). 
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We now discuss formulas (59) and (60) for a wider range of parameters. 

Remark 2 (On the values of R and I 4 ). Equations (59) and (60) can be extended to 
all physically relevant values of km (and ke) with iRa > 0 (Remark 1), by computing 
— [kpY W{kp/w, l;tr7) according to the following rules. These rules are con¬ 
sistent with the starting formulation of the Sommerfeld-type integrals (Section III). 


For zu = kmp: In regard to the residue contribution in W(0,1; kmp), set 

‘-H[^\kmP) = ± 






kmP 


kn 


upoa kn 


H\ \kmp) , (Ola) 


where ^km > 0 and the upper (lower) sign holds if the pole at km lies in the first (sec¬ 
ond) Riemann sheet; see (32). For the remaining terms in W, set a /{kmpY — {kpY/[kmP) 
a/1 — k 2 /k)n_ and choose the branch so that limfc^^oo a/1 — kS/k^n, = 1 for fixed k; thus, 

tf \eikl\<l. (61b) 

n=0^ ^/n 


• For w = kep: Set 


a /{kepy - {kpY = -i\/{kpY - {kepY = ±iujpoap /2 ; (61c) 


the upper (lower) sign holds if the pole at ke lies in the first (second) X-Riemann sheet. 


The sign in (Ola) controls the appearance or absence of the TM surface plasmon (Def¬ 
inition 1 ). If the pole at X = km does not lie on the top Riemann sheet, e.g., when k > 
and 'Isa < 0, contribution (Ola) causes cancellation of the residue from Iy,{kmP) in each rel¬ 
evant field component; then, the TM surface plasmon is absent. An analogous effect comes 
from (61c) in regard to TE polarization. Equation (61b) warrants that the primary field (in 
the absence of the layer) of the dipole is recovered as a —)■ 0 regardless of the location of the 
pole at km- See Section IVC for consequences of this solution. 


C. Field components: Evaluation 

The field components are computed by substitution of the formulas for integrals (n = 
1, 2, 3, 4) and {n = 1, 2) from Section IVB into (37)-(39) and (41)-(46), bearing in 
mind Remarks 1 and 2. To render the resulting expressions physically transparent, we 


23 



subsequently apply the condition |a;/ioO'| \k\ which warrants that \k/km\ 1 and \{ke — 

k)/ke\ 1; the requisite series expansions can be approximated by their first few terms for 
all distances p from the source. In Appendix B, we provide alternate evaluations of the field 
components for sufficiently small and large distance, p, from the source. 

Remark 3 (On the pole at km)- We assume that the parameter a is ehosen so that the 
pole at \ = km lies in the top X-Riemann sheet (see Remark 1), unless we state otherwise. 


1. Vertical dipole 


The p-component of the electric held deserves special attention, since it vanishes at z = 0 
in the absence of the layer.By using this component, we illustrate the role of the location of 
the pole at km in the Riemann surface associated with the dual (A) variable. By (37), we write 



Notice the term proportional to [2k‘^/{oopoa) — i\/k(^ — k‘^]H[^\kmP)-, which expresses the 
TM snrface plasmon; see Remark 2. By (61a), in this term = ±i2k'^/{upoa), 

depending on the location of the associated pole. On the other hand, in regard to the second 
line of (62a), k(^am ~ k^^ll — /c^/(2A;^)] ~ if [/CmP 3> |/i:p. Snppose that the pole at km 
is present in the top A-Riemann sheet, if, say, Aa > 0 and k > 0 as outlined in Section IIIC. 
By retaining two terms in the series expansion of (62a) and enforcing |a;po<^| ^ |/i:|, we obtain 


n(vd) I capoO' 


/< 


3i 


Up 


27r 


I 


X 1 + 


4 \ k'^p^ 

(cjpoo-f 


+ fTT , (kmp) + 


8F 


k^p^ k^p"^ j {ujpQaY 

1 - ^(KiikmP) - YiikmP)) + 


{ojpQCjf 


{kmPY {kmPY 


(62b) 


The TM surface plasmon contribution has survived here as the term proportional to 
Hi^\kmP)', see also Appendix B. In contrast, if km lies in the second Riemann sheet, if, for 
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example, Q'a < 0 and /c > 0, the above residue contribution disappears (see Remark 2), viz., 


(vd) (^1^0 j 


3i 


2tt 


k'^p'^ 


I + 


Ak^ 


X 




{ujpoaf 
1 3 


1 + 


8P 

, \ojpQ(j\ |A;| . (62c) 


A similar calculation can be carried out for the other held components of the vertical 
dipole. We henceforth emphasize mainly consequences of Remark 3. 

Next, we turn our attention to the magnetic held. By (38) and Remark 3, we compute 


^(vd) _ ^h'O 


P7 = -P<!2( 


ikp 


k‘^k„ 


dTT 


1 + 


2ik‘^ 


P" P 


upoa \ upoa ^kl^ - k'^ 


Hi^\kmP) 


Aik^ 


1 - ^(]HIi(pp) - Yi(pp)) 


(a;po<^)^ 

(-1)'^' p-e<' 

1=0 
ipo 


l\ 


(k'^Y (l-e<\ 

[klJ C ) 


C=kp 


-ni4 + -lP‘’' + 4>^ 

dTT 1 I p2 p 


k^ 


H[^\kmp) - i- 


Ak^ 


— YlikmP)) + 


klrP^ 


{upoaf ^ ^ (upoay 

\upoa\ < \k\ . 


1 - ^{M.i{kmP) 


(63) 


The terms with the factor in the fourth line (but not on the hrst line) of (63) correspond 
to the primary held of a ; 2 -directed dipole in free spaeef^ the remaining terms in the second 
equation amount to the scattered held due to the presence of the thin layer for \ujpQa\ |A;|. 
Notice the surviving residue contribution expressed by the term containing H^\kmp)] this 
is identihed with the TM surface plasmon via Dehnition 1. In Appendix B, we provide 
alternate calculations for in the disparate cases with \kp\ 3> 1 and \kmp\ 1- 
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The terms proportional to in the fourth line of (64) amount to the primary held of 
the dipole.^® The TM surface plasmon contribution is expressed by the term proportional 
to H‘^\kmP)- If ^ \k\ and \kmp\ = 0(1), the primary held is subdominant to this 

contribution provided \kp\ 3> \{ujpQa)/k\ (in addition to having \kp\ 1). See Appendix B 
for asymptotic evaluations of in the far- and near-held regimes. 


2. Horizontal dipole 

Next, by (41)-(46) we compute the held of the horizontal dipole. We address the ; 2 - 
components hrst, since each of these helds involves only one type of polarization, as indicated 
by the presence of only one denominator {V or Q) in their Fourier representations for ki ^ k 2 - 
By (41), involves only TM polarization and is computed as 



Notice the term (2/c^ — iupoa^/k^^ — k‘^)H[^\kmP), in conjunction with Remark 2. By 
Remark 3 and \ujpQa\ \k\ (i.e., j/CmP 3> \k\^), we keep two terms in the series expansion 
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and readily compute (see also Appendix B) 


(hd) , 

E2z -^ COS0< 


lufi^a 


2k‘^km 

(ca/iocr)^ 


1 3i 
+ 


I - An 


k^ 


k?'p^ k^p^ k^p"^ 

1 ~ 7r(EIi(^mP) ~ Yi{kmP)) + 


{(jpocr)^ 

1 3 


Hi^\kmp) 




(65b) 


In the last formula, the TM surface plasmon is evident as the term proportional to H^\kmp) ■ 
In contrast, if km lies in the second Riemann sheet (see Remark 2), then 


4“> ~ 


tupo 


cos 0< 


tupoa 


2k‘^km 

{upoaY 


1 3i 

+ 


pkp 


k^p"^ k^p^ k'^p^^ 


1 (Wi{kmP) Rl(fcmP)) + ,2 2 L4 ^4 

^mr ^mr 


(65c) 


We next compute the ; 2 -component of the magnetic held by resorting to Remark 2. Equa¬ 
tion (44) yields 


^2z 


^ sin0<| -2i — \ ^ uptake H^\kep) - 2ike\/kl - P 


ik 


in 


(1)/ 


Stt 


X 


p 


p 




z=p 


(66a) 


which involves TE polarization; and p are dehned by (56a) and (58). The term propor¬ 
tional to e^^P in the hrst line of (66a) corresponds to the held of the x-directed dipole in 
free space. Beware of the factor ^Jkl — k"^ in (66a); see also Remark 2 with w = k^p- To 
simplify formula (66a), consider \upQa\ \k\, by which ke ~ k. Thus, (66a) is reduced to 


j3(hd) *ho . 


-2i[ — - — ) e*‘'’ + ^-uiioak 

p^ pi 2 


+ 2ik^/W^^e^^P , 


H't\keP) - 2 


^kl - k^ 

upoa 


J{kep-TTlA) 


nkeP 


(66b) 


where T is dehned in (56a). If we adhere to Remark 3, we expect that (66b) should not 
manifest the contribution, //{^^(/cep), of the residue at the pole A = To demonstrate its 
elimination for large enough distances in this component, we use the asymptotic expansion^^ 
{z) ~ ^j2/{nz) which provides a reasonable approximation for H^\kep) 

for all but small values of \kep\- Recall that = —iupoa/2 if ke does not lie on the 
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top Riemann sheet of \/k‘^ — A^. Thus, for sufficiently large \kep\ ~ |/cp|, (66b) becomes 


B. 


(hd) _ *ho 


1 ik 


2z 


= sin0<(-2i + ujpoak , 


n 


Jkp 


Stt 


P" P 


kp 


(66c) 


which manifests the e*^^-wave behavior in view of (56b); see Appendix B. In contrast, by 
the hypothetical scenario when ke lies in the physical (top) A-Riemann sheet, (66b) yields 


,(hd) _ *ho 


ik 




,i{ksp—'n-l4:) 


=^sin0<j-2.^---je- + a;poaA;y-e 

in view of Remark 2. 

In regard to the p-component of the electric held, formula (42) entails 


(66d) 


pi(hd) _ 

^2p - 


(^Po 

An 


COS(f)< 


k^ 


k^\ I 


' ^ klJp^V kl kljkp^ 


Pkp 






X E (0 (P^) 


z=p 


— 2ni- 


k^ 


[upoa) 


2k^ 

(jjpoa 


-iVkl-k^ 


YiikniP)) + ^ 

rvn 


1=0 


21, l\ 


H[^^\kmp) + i 


kl 


Ak"^ km 


/ \ 0 

-2(H 



1 

V c } 

C=fcp_ 

/ 


(67a) 


where the prime denotes derivative with respect to the argument and EI)^( 2 :) = ^[]HIo( 2 :) — 
1012 ( 2 :)] + r(5/2)“^(4pr)“^2:. By imposing \ujpQa\ \k\ under Remark 3, we keep one term 
from each of the above series expansion and hnd 


E. 


(hd) 

2p 


(^Po 

47r 


cos 6 < —2\ 


1 i 
+ 


kp^ k'^p^ 


. upoa 
+ ^^~ri — 

Akp 

8k^ 

— in 


(cnpoo-)' 


hP(Kp) - 


HP'iKnP) + 


Jkp 


J{kep-3nl4) 


nkeP 
Ak^ 
{ujpoaf [ 




n{W^{kmp) - Yl{kmp)) + 


{kmpY_ 


(67b) 


The primary held of the dipole is singled out in the hrst line of formula (67b), while the 
remaining terms of this formula express the scattered wave; see Appendix B. In the scattered 



wave, there are two noteworthy contributions. The TM surface plasmon contribution is 
expressed by the term proportional to {kmp)- The other contribution, H\^\kep)i is 
cancelled out for large enough \kep\ in a fashion similar to that shown for above. 

By (43), we likewise obtain the exact formula 




sin0 


k , 

L^e 


Ak^k 


21,2 


1 ik!^ 1 


kl p\ 


Att 

+ ^-^upo<^B[^^'{kep) + \/kl- 


^ikp 


X 




3 + 4/^ {i/2Y / Ay_ k\ "" ^ \z~^/‘^Bn{zV\ 

^V2A(1-2/)(3-2/) l\ \ k^ ) dz^^^ 

r)2 Vca/ioO- J kmP 


[upoa 

oo 

yi{kmp)) + ^ 


Aik'^ ar, 
{upoaY p [ 


1 - - 
2 


1=0 


U2 



- 

\ 

V c ) 

C=kp. 

J 


\z=p 

iWiikmP) 

(68a) 


By Remark 3 and \ujpQa\ ^ \k\, this formula is simplihed to 


d“’ 


(^Po . 


1 


47r 


sin0< — + + 


ITT 

+ —upoa 


p kp"^ kP'p^ 

H^i'^\k,p) -1 


Jkp 


J{kep-3TTf4) 


nkeP 


+ 


nr , ^„ T^H[^\kmP) 


(ujpoa)^ kmP 

8k^ 1 


{ojpoaf kmp L 


l-^(ei(A;,„p)-Fi(A;^p)) + 





(68b) 


In the first line of approximate formula (68b), we display the radiation held in the absence 
of the layer. The additional terms of this formula account for the scattered held. See 
Appendix B for alternate computations based on asymptotics. 
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By (45), the p-component of the magnetic held reads 


p(hd) _ ho . 
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3 + 4/2 {i/2yfke-k\ 




^V2A (l-2/)(3-2/) 

By Remark 3 and condition \upQa\ |/i:|, the last expression becomes 


. (69a) 
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(69b) 


see Appendix B for alternate derivations of the far and near held. 

We conclude this section with the computation of the 0-component of the magnetic held. 
By (46), we obtain 
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(70a) 
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Suppose that \uJiJ,Qa\ By Remark 3, the preceding formula is reduced to 


jj(hd) ho , . , * , 

^26 ^ T- cos — + 


^ikp 


An 


kp^ k‘^p^ 




{ujpoaY 

in {up^aY 
8 keP 


{upoay 


^{m[{kmp) - Ylikmp)) + 


i7«(M - 


p{kp-‘i-Kl4) 


nkeP 




(70b) 


In Appendix B, we compute this held component by asymptotic methods, in the near- 
i\kmP\ ^ 1) and far-held {\kp\ 3> 1) regimes for \ujpQa\ \k\] the results are in agreement 
with (70b). 


V. DISCUSSION 

In this section, we discuss the derived formulas of Section IV under Remark 3 in the case 
with \kmp\ 3> 1 and kp = 0(1), by assuming that \ujpQa\ \k\. An issue that we address 
is the relative strength of the TM surface plasmon contribution to the electromagnetic 
held. A practical consideration that motivates this study is that layers of two-dimensional 
materials snch as graphene may sustain a TM snrface plasmon wave nnmber, km-, with 
0 < ^km 'Alkm and ikkm. ~ \km\ 3> \k\ at a certain accessible freqnency range;^'^ thus, the 
wavelength of the generated snrface plasmon can be much smaller than the wavelength of 
the incident electromagnetic held in free space. It is of interest to investigate how the TM 
snrface plasmon can prevail over other held contribntions on the layer at distances of the 
order of the free-space wavelength. 


A. Field of vertical dipole, \kmp\ 3> 1 and kp = 0(1) 


By (62b), the p-component of the electric held reduces to 


(vd) _ f _3*_ Ak^ i(k^p-n/A) 

Y 27r ^ 4 \jc^p‘^ kP'p^ k'^p'^J (ujpQa)'^\j kmP 


(71) 
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if \kmp\ 3> 1, by use of asymptotic expansions (54) and v/ 27 Me*(^-^^/ 2 -^/ 4 ) as 

z ^ oo with u = 1. Hence, the TM surface plasmon dominates in this regime if 

l^h'OO'l ^-{2/7)‘ikm/\k\ 


k 


< e 


In the same vein, by (63) we compute the 0-component of the magnetic field: 
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^ i{kmp+Tr/ 47 ) 
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dTT I kpj^ {upoay\l 2kmP^ j 

Evidently, the TM surface plasmon contribution dominates in this component if 

^^ 0^1 ^-(2/3)Sik^/lkl 
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Equation (64) for the ; 2 -component of the electric held yields 
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which implies that the TM surface plasmon is dominant if 
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< e 


-{2/5)7sk^/\k\ 


( 72 ) 


( 73 ) 


( 74 ) 


f 


( 75 ) 


( 76 ) 


Thus, for kp = 0{1) the TM surface plasmon dominates in all held components of the 
vertical dipole if (74) holds. 


B. Field of horizontal dipole, \kmp\ 3> 1 and kp = 0(1) 

First, we address the ; 2 -components. By (65b), the ; 2 -component of the electric held reads 

ik"^ 

+ — - — 1 ^ ^ J. . (77) 


i-,(hd) ^po , j Ujp^a 
^2, ~ — COS0' 


271 4 \ /c2p2 ' /c3p3 /^4p4 j {uJPqCtP \ k^P 

For kp = 0(1), the pole contribution dominates under condition (72). In contrast, by (66) 
the ; 2 -component of the magnetic held does not carry any TM surface plasmon contribution; 
thus, we omit its evaluation from this discussion. 

The remaining components are more richly structured. By (67b), we have 
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(78) 
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If kp = 0{1), all terms in the second line are relatively small. The residue contribution of 
the pole at km dominates over the primary held if (76) holds. 

Equation (68b) for the ^-component of the electric held gives 





gifcp _ 


8k^ 


(upoa)^ kmp V kmp 


271 


Pikmp-Trli) 


i{kep—STrl4) 




(79) 


which is dominated by the TM surface plasmon at distances with kp = 0{1) provided (74) 
holds. 

We now focus on the components of the magnetic held. By (69b), we compute 
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All terms of the second line are negligible; in contrast, the TM surface plasmon can dominate 
over the e*^^-wave of the hrst line under (74). Similarly, formula (70b) is reduced to 


jj(hd) Pok , 
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1 i 
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(81) 


which suggests the dominance of the TM surface plasmon if (76) holds. 

Therefore, on the basis of the approximate formulas of Sections VA and VB, we reach 
the following conclusion. 


Proposition 1 (On dominance of TM surface plasmon contribution). Suppose that the 
pole at X = km is present in the top X-Riemann sheet; and \upoa\ <A \k\ which amounts to 
C •= Wl^/JpJWl ^ 1? where k = Uy/poi. Then, the TM surface plasmon is dominant at 
distances p of the order of the wavelength in free space if condition (74) holds. In particular, 
if k > 0 (e > 0) and arga = 7;/2 — 6, 0 < 5 < 7r/2, condition (74) reads 


( <A exp 



(82) 
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VI. CONCLUSION 


In this paper, we computed by analytical means the electromagnetic field of a vertical 
and a horizontal dipole located near a thin film with scalar, complex conductivity, a, inside 
an unbounded medium of wave number k. This setting provides a minimal yet nontrivial 
model of wave propagation on a two-dimensional material such as a graphene sheet. 

Our starting point was the Fourier-Bessel representation, or Sommerfeld-type integrals, 
for all field components. In the particular case when both the dipole and the observation 
point lie on the layer, we showed that the fields can be expressed in terms of geometrically 
convergent series and known transcendental functions such as the Fresnel integrals, and the 
Bessel and Struve functions. We simplified this result considerably when |a;/ioO'| ^ \k\ and 
derived analytic formulas that hold for practically all distances, p, from the source. These 
formulas connect smoothly the near field with the far field of the dipole, as we verified via 
independent asymptotic evaluation of the requisite integral representations. Our treatment 
provides an analytical description of a particular case of the dyadic Green function associated 
with the geometry of the conducting film; and can expand insights obtained from previous 
numerical or semi-analytical approaches for related problems. 

Depending on the frequency u and the parameters a and k, the electromagnetic field of 
the dipole may manifest a contribution from a pole analogous to the TM surface plasmon 
of plane waves.^ We demonstrated explicitly how this contribution may be manifest or 
suppressed in each relevant component of the electromagnetic field for all distances from the 
dipole source. 

The present work admits several extensions and also points to pending issues. For ex¬ 
ample, we have not computed the field components at points away from the film; this task 
would require an elaborate asymptotic evaluation of the requisite Fourier-Bessel integrals 
by imposing some restriction on \z\. The possible anisotropy of the thin-film material, e.g., 
black phosphorus,^® where a is a tensor, deserves attention; this anisotropy dramatically 
affects the dispersion relations for surface plasmons and the analytic structure of the elec¬ 
tromagnetic field produced by current-carrying sources. Another problem is the excitation 
of surface plasmons by a resonant receiving antenna placed on the material sheet,which 
requires solving an integral equation for the associated current distribution; this task is left 
for future work. 
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Appendix A: On the normal field components 


In this appendix, we discuss the conditions on the field components normal to the plane, 
z = 0, of the him. Such boundary conditions were not explicitly invoked in the derivation 
of the Fourier-Bessel integral representations (Section III). 

First, we describe the (effective) surface charge density of the him which produces a jump 
of the normal component, ijEjz, of the electric displacement held = ejEj. By (4b), the 
thin layer bears the surface current density J* = aEji = a{Ej — (Ej ■ 62 ) 0 , 2 } a.t z = 0, which 
by the continuity equation is associated with the surface charge density 


Qs{.x,y) = -{i/uj)S/s ■ Js ^^ 0 = -{icr/u)Vs ■ Ey ; V,, = V - (V ■ 62)62 . 


By Gauss’ law in region j, Vs-Ey = —dEj^/dz for {x, y, z) 7 ^ (0, 0, a); allowing 2 : to approach 
0 from region j, we conclude that the normal derivative, dEjz/dz, of Ejz is continuous across 
z = 0. Thus, 


0s = 


ia dEiz 


u dz 


ia dE 2 z 


2=0+ 


u dz 


(Al) 


2=0- 


where writing z = 0^ (;2 = 0 “) implies that = 0 is approached from above (below). 

The boundary conditions for the normal held components consist of: (i) the jump of 
the 2 :-component of the electric displacement held due to Qs] and (ii) the continuity of the 
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2 ;-component of the magnetic field. Accordingly, we write 


[k^Ej^] : = {klEu - klE2z)\z=o = Qs 


dE., 

lu/ioa 


[Bjz\ — — B2 z)\z=0 — 0 . 


(A2a) 

(A2b) 


Next, we proceed to demonstrate that (A2) are satished by the integral representations 
of Section III for an elevated dipole (a > 0), as expected. 


1. Vertical dipole 


By differentiation in of (14b) for 0 < 2 : < a and (16b) for < 0, we obtain 


dE,z 


dz 


j=0+ 


^fdAA3Ah,. 


OE. 


2z 


dz 


(A3) 


z=0- 


which shows the continuity of the normal derivative of Ejz- On the other hand, we compute 




] = r T J„(Ap) (1 + 

/3i y /c2/3i + ki(32 + i^/^o<^/3i/32 j 


An 


27r 


dA A^ Jo(Ap) 


kl 


k^Pi + kl(32 + i^ho<^/3i/32 


3*/3ia 


= (uji^oct) dA A^ Jo(Ap) ^ , 

which, by (A3), is equal to iufXQa{dEjz/dz) at z = 0. Thus, condition (A2a) is satished. 
Note that (A2b) is trivially satished in this case because Bj^ = 0. 

2. Horizontal dipole 

By diherentiating (24c) and (26c), we verify that dEjz/dz is continuous across z = 0: 


dEu 


dz 


dE- 


’:= 0 + 


/'°°^\ \2 T M A ^1^2 i/3^, 

■ —cos^/ dAA J.(Ap)—a. 


22 


(A4) 


2=0- 


On the other hand, the jump of k‘^jEjz is equal to 

r,2p. 1 Tf\ ^ (k‘iP2-kjPi-UHo(TPiP2 ^ 1 

^ ^ ' 47r ^ kjP2 + klPi + uiioa^i^2 


H- — cos(j) / dAA^Ji(Ap) 

zn In 


kldi 
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which, after elementary algebra, is equal to iui^oa{dEjz/dz) at z = 0] cf. (A4). 

We now turn our attention to By (25c), as 0 (with 0 < 2 : < a) we obtain 


Blz\z=0+ 


11^0 

47r 


sin 



dAA^ Ji(Ap) 


1 

Ji 



/32 — Pi + cofiQa 
P2 E Pi UflQ(J 


g*/3ia 


which is equal to B 2 Z at 2 : = 0 by (27c). This shows the expected validity of (A2b). 


Appendix B: Asymptotic evaluation of Sommerfeld-type integrals 


In this appendix, we derive the near- and far-held of each dipole by asymptotic evaluation 
of the Sommerfeld-type integrals under Icapo^'l <^ |/^|- The results of these calculations are 
compared to the exact formulas of Section IV C. We assume that the pole at \ = km is 
present in the physical Riemann sheet throughout this appendix (see Remarks 1-3). 

In the following analysis, we will need the relations^^ 

Mz) + J 2 {z) = ^ Ji{z) , Jo{z) - J 2 {z) = 2J[{z) , (Bl) 

along with the leading-order asymptotic formula^^ 



Juiz) ^ \ — cos (;2 — z/7r/2 — 7r/4) as | 2 :| —)■ cx) , |arg2:|<7r 


TTZ 


(B2) 


In addition, we will make use of the following result from Ref. 26. 
Lemma B.l (On the Fresnel integrals). The integral 

poo ^i{ke-k)pT 


3 := 


dr 


TT 


, (B3) 


^/T{uJHoa + i2V^y/ke — k^/r) '/ky/ke — k 

where yV > 0 for r > 0, B{z) is defined in (56a), and po = —{ojpiQafip/ifik) [see (36) 
and (58)y. 

Proof. By Ref. 26, hrst write 


3 = 


dr 


^i{ke-k)pT 


^00 

+ / dr 


^-i{ke-k)pT 


lo y/T{ufj,oa + i2y/^y/ke — k^/T) Jo —i^/T{ufJ,oa + 2y/2k^ke — ky/r) 

Then, apply r h-)■ <^ with y/g = 2iy/^y/ke — k{ujpLo^)~^y/P fo hnd 
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where 


poo ^iz<; poo 

3„(z) = -j^ dc ^(^^1)+/^ * Vf(Vf+i) 


This 3o{z) satisfies 


^ - i3o{z) = -2e , lim ^ 0 ( 2 :) = 0 , 

dz y z 2 ^-+oo 


(B4) 


with solution 3o{z) = 2Ti\/2e which concludes the proof. □ 

1. Far field: \kp\ 3> 1 with |a;/ioO'| \k\ 

In this regime, the major contribution to integration in the Sommerfeld-type integrals 
comes from the singularities at A = /c and \ = km- The Bessel functions in all integrands 
are expressed in terms of -Jui^Xp) for z/ = 0, 1 by use of (Bl); Ju{Xp) is in turn replaced by 
its large-argument approximation according to (B2). 

a. Field of vertical dipole on film 


Taking the limit as a j, 0 with 0 < < a of (14a) and then invoking (B2) readily yield 


Eip = --^l d\\\h{\p) 
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47t ‘. V 2k'^ + Uipou^/k'^ - 

In this integral, we single out the contributions from the simple pole at X = km and the 
branch point at X = k, viz.. 


(vd) lojpo 1 2 

^ip -27r*/c^ 


dTT 




I 


_ ^^{^rnp 37r/4) 

irk^fYlP A— 


2/c^ — up^a^/kX — A^ 
(a;/ioO-)2(2A) 


-|- ik^ / dr {1 + irf 

Jo 
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X 


7rkp{l + ir] 


^Hkp-O-KjA) ^-kpT 


where for the last integral we appropriately deformed the integration path in the upper 
A-half-plane, set A = k{l + ir) and integrated on each side of the positive real r-axis. The 
major contribution to integration comes from the vicinity of r = 0 with width 0{{kp)~^)-, 
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the corresponding integrand yields a term of the order of ufj,Qa/k. By enforcing the condition 
Ica/ioO'l \k\, we reduce the preceding formula for to 
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(B6) 


kmP ' ‘2k k^p"^ 

V 

in agreement with (62b) which comes from the exact series expansion for this component. 
In the same vein, by (14b) and asymptotic formula (B2) for z/ = 0 we have 


p(vd) _ 
^Iz ~ 


upo 


2nk^ Jo 

|fcp|>l ojpo 


dA Jo(Ap) 


A^ /c^ + upoay/k'^ — A^ 


Ank"^ 


dA 


— A^ 2A;2 + ojpoay/k'^ — A^ 
+ ojpoay/k'^ — A^ 


A3 


„i(Ap-7r/4) 


a/A:^ — A^ 2 A;2 + ijjpQa\/k‘^ — A^ V "^Ap ^ ^ 

By singling out the contributions from two singular points, at A = A;^ and \ = k, and 
additionally imposing \upo<^\ ^ |A^|, after some algebra we hnd 
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by expanding the integrand at r = 0. The preceding result is in agreement with (64) of the 
exact solution for Icapo^'l ^ \k\ and |A;p| 3> 1. By a similar computation, (15) yields 
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which, alternatively, results from (63) of the exact solution if |a;/xoO'| ^1^1 and \kp\ 3> 1. 


h. Field of horizontal dipole on film 


In this case, the contribution from the branch point at X = k should be affected by the 
proximity of the pole at \ = ke if |(/Ce — k)p\ < 0{1). This effect is taken into account via 
the integral of Lemma B.l, where {ke — k)p ~ po enters as a parameter. This situation is 
typical in radiowave propagation on the boundary separating media of signihcantly different 
wavenumbers.^® 

In view of (Bl) and (B2), the limit as a 0 and 2 : t 0 of (26a) yields 
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There is only one residue contribution, component, which comes from the 


last integral and concerns the simple pole at \ = km'- 
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The contribution from the branch point at A = /c is singled out as 
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For the hrst integral above, we have set X = k + {ke — k)T'\ whereas for the second integral 
we deformed the path in the upper A-half-plane and then set A = A;(l + ir) as before. In each 
of the transformed integrals, the major contribution arises from a small neighborhood of the 
origin. The value of the hrst integral crucially depends on which term in its denominator 
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prevails when {kg — k)pT' = 0(1); by this interplay, the parameter po emerges naturally. In 
contrast, in the second integral such interplay is not present; one may simply Taylor-expand 
the integrand in powers of ojpou/k and keep the leading-order term. To subtract the primary 
field contribution in the hrst integral, we split each term, evaluated near r' = 0, as 


ujp,Q(j 


2i\/^y/ke — /cV^ -|- oopoa 2i^/^\/ke — ky/^ {2i^/^^/k^^^\/~T'+ ufioo) 

and carry out the integration without further approximations for these terms. Thus, we hnd 
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where 3 is given in Lemma B.l. Thus, by rsj ^£11 + we compute 
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which is in agreement with (67b) for \kp\ 3> 1. 

The remaining components can be approximately computed by the same methodology, 
although the algebraic details are different in each case. By (26b), we have 
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By separating the pole and branch-point contributions, we hnd 
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which, by Lemma B.l, is in agreement with (68b) if \kp\ 3> 1. We note in passing that 
the e*^^/(A;^p^) term in the hrst line of (B13) arises from the (TM-polarization) integral 
containing {2k‘^ +upoa\/k‘^ — A^)“^, is small compared to e*^^/p in this far-held regime and, 
thus, should be dropped to leading order in kp. 
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By (26c), the ; 2 -component of the electric held does not contain any TE polarization term 
and thus lacks the Fresnel integrals. By similar asymptotics, we reach the following result. 
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cf. (65) from the exact solution, recalling that km ~ i2k’^/{ojpoa). 

Next, we turn attention to the components of the magnetic held. By (27a), we have 
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in the third line we replaced — A^/ {2\/kA — X^ + up^a) by 1/2 — {up^a/2) / {2\/kA — A^d- 
upoa) and integrated out the constant. For the hrst integral, we deform the path in the 
upper A-half-plane, picking up the residue from A = km, and then evaluate the contribution 
from the vicinity of \ = k. For the second integral, we set A h-)■ r with r = {\ — k)/{ke — k) 
and evaluate the contribution from the vicinity of r = 0, ending up with Fresnel integrals 
via Lemma B.l. Accordingly, we derive the expression 




'2p 


Att 


. { Airik"^ 1 

I {ojpoaY kmP V Tr/cmP 


gi(fc™p-37r/4) li^PoCr 


k'^p^ 


pkp 


A \ kp 


(B16) 


the e^^^/{k‘^p^) term in the parenthesis of the hrst line comes from the TM-polarization 
integral and is neglected. The ensuing asymptotic result can also be obtained from (69b) in 
the far-held regime, \kp\ 3> 1. 
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In this vein, by (27b) we compute the following expression [cf. (70b)]. 
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Here, the branch-point contribution for each polarization (TE and TM) integral carries a 
free-space term proportional to [a;poO'/(A;p^)]e*^^, besides the Fresnel integral contribution. 

We conclude this subsection with the approximate computation of the ; 2 -component of the 
magnetic held by similar means. This component expresses only TE polarization and, thus, 
lacks any pole contribution (by Remark 3). By (27c) and similar algebraic manipulations 
as above, we hnd 
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in agreement with (66c). 


2. Near (static) field: \kmp\ 1 with \ujpQa\ \k\ 

In this case, the major contribution in each Sommerfeld-type integral comes from A = 
C>(l/p) 3> \km\ 3> \k\. Thus, we make the approximations 

a//c^ — A^ ~ f A , 2/c^ -|- upoax/k'^ — A^ ~ iXupocr , 2a/— A^ -|- upoa ~ 2iX , 

where A > 0; and exactly evaluate the ensuing integrals containing J^{Xp). 
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a. Vertical dipole 


In the limit a 0 with Q < z < a, (14a) yields 
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Ank^upQa Jq r/ 2 ^^ oj^qO ’ 

which agrees with (62b) in the regime where \kp\ |/CmP| ^ 1 and \upQa\ \k\. In this 
regime, the factor as well as terms containing the Struve and Neumann functions must 
be expanded out in powers of kp and kmpf^ after several mutual cancellations, it is found 
that the term {kmp)~‘^ in the second line of (62b) provides the leading-order contribution. 
Likewise, (14b) furnishes the integral 

(vd) iujpo wx X iujpol ^ r\ 7 MM iujpo I , . 

via integration by parts; this result is in agreement with (64) where the leading-order con¬ 
tribution comes from the last term. The remaining component is provided by (15), which 
is reduced to the following expression [cf. (63)]. 
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h. Horizontal dipole 

In the static held of the horizontal dipole, the two types of polarization may provide 
contributions of different orders of magnitude. By (26a), we have 
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with dominance of the term pertaining to TM polarization. This result agrees with (67b), 
in which [kmp) and Yl{kmP) provide the leading-order contribution after cancellation 
of the 0{1/p^) terms. By (26b), the 0-component of the electric held is 
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in agreement with (68b) where the leading-order term comes from the expansion of 
H^\kmP) and {kmP)~^yi{kmP) for \kmp\ 1 under \ojpQa\ \k\. The remaining 
component of the electric held is evaluated on the basis of (26c), viz., 
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271 (jjpoa p^ 

This result smoothly connects with (65b) where the leading-order contribution comes from 
the {kmp)~‘^ term of the second line. 

We proceed to the computation of the magnetic held. By (27a), we write 
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This formula also results from (69b); the leading-order term comes from the expansion of 
{kmP)~^H[d(kmP) and {kmP)~^Yi{kmP)■ By (27b), we obtain 
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All p2 

Now compare this approximation to (70b); the latter reduces to the above expression by 
expansion in powers of kp and kmP- The remaining component of the magnetic held is 
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This formula manifestly agrees with (66b) if |A;p| 1 and \ujpQa\ \k\. 
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